Decoherence in Adiabatic Quantum Computation 
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We have studied the decoherence properties of adiabatic quantum computation (AQC) in the 
presence of in general non-Markovian, e.g., low-frequency, noise. The developed description of 
the incoherent Landau-Zener transitions shows that the global AQC maintains its properties even 
for decoherence larger than the minimum gap at the anticrossing of the two lowest energy levels. 
The more efficient local AQC, however, does not improve scaling of the computation time with 
the number of qubits n as in the decoherence-free case. The scaling improvement requires phase 
coherence throughout the computation, limiting the computation time and the problem size n. 



The adiabatic ground-state scheme of quantum com- 
putation [j], 0] represents an important alternative to 
the gate-model approach. In adiabatic quantum com- 
putation (AQC) the Hamiltonian Hs of the qubit reg- 
ister and its wave function \ip) undergo adiabatic evo- 
lution in such a way that, while the transformations of 
\ip) represent some meaningful computation, this state 
also remains close to the instantaneous ground state 
\ipo) of Hs throughout the process. This is achieved 
by starting the evolution from a sufficiently simple ini- 
tial Hamiltonian Hi, the ground state of which can be 
reached directly (e.g., by energy relaxation), and evolv- 
ing into a final Hamiltonian Hf, whose ground state pro- 
vides the solution to some complex computation problem: 
Hs—[1 — s(t)]Hi + s(t)Hf, where s(t) changes from to 
1 between some initial (tj=0) and final (tf) times. 

The advantage of performing a computation this way, 
besides its insensitivity to gate errors, is that the energy 
gap between the ground and excited states of the Hamil- 
tonian Hs ensures some measure of protection against 
decoherence. This protection, as partly demonstrated in 
this work, is not absolute. Nevertheless, it allows for the 
ground state to maintain its coherence properties in time 
far beyond what would be the single-qubit decoherence 
time in the absence of the ground-state protection. This 
feature of the AQC remains intact [|| even if the deco- 
herence strength and/or temperature is much larger than 
the minimum gap. 

In general, the performance of an adiabatic algorithm 
depends on the structure of the energy spectrum of its 
Hamiltonian Hs- Here we consider a situation, which 
is typical for complex search and optimization problems 
Q, when the performance is limited by the anticrossing 
of the two lowest energy states. The minimum gap g m 
between those states shrinks with an increasing number 
n of qubits in the algorithm, although the exact scaling 
relation is not known in general. In an isolated system 
with no decoherence, the limitation is due to the usual 
Landau-Zener tunneling at the anticrossing, which drives 
the system out of the ground state with the probability 
given by the "adiabatic theorem" . Different formulations 
of the theorem all give the computation time as some 
power of the minimum gap: tf oc g"" 5 [1, [f|. 

The main assumption behind the adiabatic theorem is 
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FIG. 1: Broadening of the energy levels of a closed system (a) 
due to coupling to an environment made of (b) a single two- 
state system, or (c) infinitely many degrees of freedom with 
a continuous energy spectrum. In general, the coupling to an 
environment splits a single anticrossing into M 2 anticrossings 
where M is the number of environment energy eigenstates. 
For environment with a continuous spectrum the anticrossing 
turns into a continuous transition region of width W . 



that there exists a well-defined energy gap between the 
two lowest energy states of the system. In a more realistic 
case with decoherence, however, the energy levels of the 
qubit register are broadened by the coupling to environ- 
ment, as illustrated in Fig.[T] Even the simplest environ- 
ment, e.g., a two-state system, splits a single anticrossing 
of the two qubit levels into four anticrossings with smaller 
gaps (Fig.[T]3). An environment with a continuous spec- 
trum turns the anticrossing point into a continuous re- 
gion of some width W (Fig. [TJ:) within which incoherent 
tunneling between the two qubit states can take place. 
Thus for such typical models of environment, the gap no 
longer exists in the "qubits+environment" system. The 
broadening W is directly related to the decoherence time 
of the qubit states. Any uncertainty W in the energy 
of an energy eigenstate makes the accumulated phase of 
this state also uncertain in time Td CC oh ~ 1/W. Since 
the broadening W typically increases with the number 
of qubits, while the minimum gap g m decreases, the re- 
alistic large-scale system will eventually fall in the in- 
coherent regime W ^> g m . This means that studies of 
the adiabatic theorem do not apply to such realistic sit- 
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uations and therefore new ways of understanding AQC 
performance become necessary. One possible approach 
towards this goal is to generalize the adiabatic theorem 
to open quantum systems @. 

In this paper, however, we study the evolution of an 
adiabatic quantum computer in the "incoherent" regime 
by developing a corresponding description of Landau- 
Zener transitions for W 3> g m - We use the model of de- 
coherence appropriate for solid-state circuits, where the 
AQC approach is particularly promising. One charac- 
teristic feature of such a model is that it should allow 
for non-Markovian, in particular low-frequency, environ- 
mental noise. Previous studies have mainly considered 
Markovian environments [1 0, B [§]• A correct descrip- 
tion of the interaction with a low- frequency environment, 
which has the strongest effect on the AQC algorithms 
[Irj| . requires a non-perturbative or strong-coupling the- 
ory of the environment-qubit interaction. 

Another feature of our "solid-state" approach is the 
assumption that the environment responsible for deco- 
herence is in equilibrium at some temperature T, and 
is sufficiently large to enforce (on some time scale) the 
equilibration among the qubit states at the same temper- 
ature. Even the low-frequency noise that dominates the 
decoherence of the solid-state qubits (see t e.g., [ill Il2j |) 
comes usually from equilibrium sources [131 ] - Previous 
studies of the AQC decoherence used models that do not 
account directly for such equilibration [f| 0, H, [13, EH • 
Since environment temperature can not be reduced indef- 
initely, for a sufficiently large system, T will inevitably 
be larger than the minimum gap g m . This means that 
in contrast to closed systems, Landau-Zener transitions 
in the presence of decoherence are intrinsically linked to 
thermal excitations out of the ground state, making it 
necessary to consider the two types of the transitions si- 
multaneously. 

Quantitatively, we introduce the decoherence as usual 
by adding the bath Hb and the interaction Hamilto- 
nian H lvA to the Hamiltonian H$ of the qubit register: 
Htota\=Hs + Hb + -ffint- As discussed above, we use the 
two-state approximation near the anticrossing, assuming 
that g m is much smaller than the energy gaps separating 
the first two from the other levels 1 231: 



-(ea z + g m a x )/2, 



(1) 



where cr's are the Pauli matrices, Q is an operator of the 
environmental noise, e = E(s—s m ) with E 3> g m defin- 
ing the energy scale which characterizes the anticrossing 
at s = s m . Independent couplings of individual qubits to 
their environments produce only the uncoupled noise in 
the two-state model (H} Q. We assume that the noise is 
Gaussian so that we do not need to specify Hb explic- 
itly. Then, all averages can be expressed in terms of the 
spectral density: 



S(u) 



dte iuJt (Q(t)Q(0)), 



Gaussian noise is expected if the environment consists of 
a large number of degrees of freedom all weekly coupled 
to the system [jfjj . 

In the regime of incoherent Landau-Zener transitions 
considered in this work, both the environment-induced 
broadening W of the two basis states of the Hamiltonian 
(fT|) and temperature T are taken to be much larger than 
g m . This means that the time (~ 1/W), during which 
the two states lose their relative phase coherence, is much 
smaller than the typical tunneling time (~ l/ffm) which 
implies that the tunneling between these states will be 
incoherent. In particular, the off-diagonal elements of the 
density matrix p of the system (fTJ vanish within the time 
Tdecoh ~ 1/W so that p reduces to diagonal elements, i.e. 
to p z = po — pi , which is governed by the usual kinetic 
equation 



-T(p z - poo), 



(2) 



where r = r i + r 10 and p^ = [T w - r i]/r. Here 
we use the standard notations: |0) and |1) are the two 
eigenstates of a z with eigenvalues =Fl, respectively, pj is 
the occupation probability of state \j), and is the rate 
of tunneling from state \i) to \j). 

The physics behind such an incoherent tunneling is the 
same as for macroscopic resonant tunneling (MKT) of 
flux in superconducting flux qubits which has been stud- 
ied experimentally [l3l| and theoretically [13] • In partic- 
ular, the transition rates have the structure of resonant 
peaks of width W in the vicinity of the anticrossing point. 
These rates can be explicitly calculated by a perturbation 
expansion in g m and assuming Gaussian noise [171 ]: 



r i(e) = 



<J7, 



dte let exp 



du e-*"*- 



(3) 



The rate of the backward tunneling is determined by the 
relation rio(e) = roi(— e). In the case of white noise, 
S{uj) — 5(0), Eq. ^ gives the tunneling peak in the 
form of a Lorentzian line-shape: 



r i(e) 



1 9 2 m W 
2e 2 + W 2 ' 



W=-S(0). (4) 



On the other hand, in the situation characteristic for 
practical solid-state qubits when the noise is dominated 
by the low- frequency components, Eq. ([3]) reduces to a 
shifted Gaussian (l7j : 



2W 2 
f du S(u) 
p ~ ' J 2n w ■ 



(5) 



For environment in thermal equilibrium, the width W 
and the position e p of the Gaussian are related by [lTj : 



W z = 2Te 



(6) 



where (...) denotes averaging over the environment. 



These theoretical results have been experimentally con- 
firmed in flux qubits [l3j]. 



3 



Let us first study the kinetic equation ((2]) in two ex- 
treme cases. In the small-T regime p^ ~ sgn e which 
implies, with the initial condition p z (0) = 1, that the 
right hand side of © is nonzero only for e > 0. This 
leads to the ground state probability 

VG = l-e-^f, (7) 
1 f (1 de 1 f°° , s de 

i = T r( e )- = -/ r 01 (c)-. (8) 

*/ Jo e */ J-oo e 

We shall see later that under relatively general conditions 
e oc 1/tf and therefore 7 is independent of tf. These 
equations assume that the range of e is large enough to 
effectively cover the whole peak of Toi, therefore justify- 
ing infinite integration limits. In particular, the range of 
e should be larger than (among other energies) the cutoff 
energy of the environment excitations. In the opposite 
large-T regime, one has |e| -C T and hence poo = in 
Eq. © for energy e within some relevant interval around 
the anticrossing point e = (this condition is made more 
precise below). The ground stare probability is then 

PG = \(l-e-^'). (9) 

Because of the thermal excitations, pc approaches 1/2 in 
the slow-evolution limit. For the intermediate T regime, 
Pg always falls between (J7J and |9j , therefore these equa- 
tions give, respectively, upper and lower bounds for the 
probability of success (see Fig. [5] and discussion below). 

An important feature of §S§ is that for uniform evolu- 
tion, i.e., e = const = u, it gives ytf — - Toi(e)de = 
wg^ n /2i', independently of S(lo), leading in the small-T 
regime to the same Landau-Zener probability ([7]) as in 
the decoherence-free case. This result extends the recent 
proofs [TH, [n| [2(| that at T = Landau-Zener probabil- 
ity is unaffected by decoherence. The physical reason for 
this is that the decoherence changes only the profile of the 
transition region while keeping the total transition prob- 
ability the same. Therefore, in the two extreme regimes, 
the ground state probabilities and © are completely 
independent of the form of the noise spectrum S(u>). 

At intermediate temperatures, on the other hand, the 
quantitative ^/-dependence of the probability pq is sen- 
sitive to the specific form of S(u>) and therefore to the 
tunneling rates. For Gaussian rates ([5]) and uniform evo- 
lution, pa calculated from Eq. © is shown in Fig. [2] 
The curves characterize the transition between the low- 
and high- ([9]) temperature limits. At small evolution 
times when tf <C 7 -1 all curves coincide, with pq = jtf 
in the linear approximation, independently of temper- 
ature T. The temperature dependence of pq appears 
only in the second-order terms in jtf. For slow evolu- 
tion, tf > 7 _1 , pg varies from 1 to 1/2 with tempera- 
ture T - see inset in Fig. [5] If the evolution is infinitely 
slow, the occupation probabilities of the states |0) and 
|1) should always reach the local thermal equilibrium. 
This, however, is not the relevant regime for the present 
discussion. In the relevant case, the rate v is compa- 
rable to the maximum tunneling rates T, and therefore 
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FIG. 2: The occupation probability pa of the ground state as 
a function of the dimensionless evolution time yt / for different 
temperatures T in the case of the Gaussian tunneling rates 
@. The inset shows the dependence of pg on T/W for jtf = 
1; 1.5; 2; 3; 5 from lower to upper curves respectively. 

becomes much larger than the tunneling rates as the sys- 
tem moves away from the resonance, so that the local 
equilibrium is not maintained. This means that, strictly 
speaking, the large-T result © is valid for any t f only for 
T 3> E. Asymptotic analysis of the evolution equation 
for the case of the Gaussian rates (JSJl shows that in the 
more interesting regime when T 3> W but T <C E, the 
ground-state probability is: 

PG = \ + ^f [lnitf]1 ' 2 - (10) 

This equation describes the increase of pg towards the 
local equilibrium at sufficiently large evolution time tf, 
and corresponds to the large-T part of the two curves 
with larger -ft f in the inset in Fig. [2] 

We now use the results presented above to discuss the 
performance of AQC in the incoherent regime g m <C 
W,T. For this, one needs to distinguish global and local 
adiabatic evolutions. In the global scheme, the adiabatic 
evolution is uniform, e = const = E/tf, and Eqs. ([7]) 
and © show that the required computational time tf ~ 
7 _1 = 2E/irg'j tl coincides with the decoherence-free case 
independently of decoherence and temperature T. Even 
if the large T reduces pc to ~ 1/2, to find correct solu- 
tion, one only needs to repeat the computation process 
on average two times. 

Global adiabatic evolution, however, does not yield the 
optimal performance in coherent AQC. Indeed, for the 
case of adiabatic Grover search [l4[ , the global adiabatic 
scheme yields the complexity of the classical exhaustive 
search, i.e., tf — O(N), where N (= 2") is the size of 
data base. In the more efficient local scheme [3], one 
takes e(t) = ag(t) 2 , so that the adiabatic condition is 
satisfied uniformly (the system slows down in the region 
of small gap) and the computation time is tf = n/ag m 
which for the case of adiabatic Grover search yields the 
optimal 0{^/N) performance. The local evolution plays 
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crucial role for the scaling analysis of the AQC [IJ, [2l|, 
[22J, although in some cases it is only assumed implicitly. 
In general, however, finding the gap g(s) is as hard as 
solving the original problem, and only in some cases, e.g, 
the adiabatic Grover search, g(s) is independent of the 
final solution and can be found a priori analytically. 

The enhanced performance of the local scheme comes 
at a price of its stronger sensitivity to decoherence. A 
qualitative reason for sensitivity of local AQC is that 
although decoherence does not change the total integral 
transition probability, it distributes it over a much larger 
energy interval W ^ g m , making it necessary to slow 
down the evolution for a longer period of time. If one 
uses the same e(t) as in the decoherence- free case, the 
average tunneling rate (|8|) is dominated by the vicinity 
of the point e = 0. Quantitatively, e = ag 2 and tf = 
Tr/affm yield (i/e) -1 = g m /irg 2 w S(e), which together 
with and ^) give 7 w r i(0) oc g 2 m . Therefore the 
computation time is tt ~ 7 -1 oc g~ 2 , which is similar 
to the performance of the global scheme with the only 
possible enhancement compared to the global case being 
a prefactor. In the case of white noise, Eq. (j4|) leads to 
7 = gf n /2W, while for the low-frequency noise, Eq. ^ 
gives 7 = y / Trj8(gl l /W)e~ w2/8T2 . Notice that in the 
latter case, lowering T with constant width W [HI does 
not shorten the computation time. 

To summarize, we have studied the decoherence effects 
on AQC due to general non-Markovian environments in 
the strong decoherence regime in which the broadening 



of the energy levels completely smears out the anticross- 
ing region. Our strong coupling treatment shows that 
global AQC remains unaffected by strong decoherence 
W > g m and is independent of the type of noise, while 
the local AQC provides only a prefactor improvement 
of the algorithm running time in this regime and does 
not change the scaling of this time with g m as com- 
pared to the case without decoherence. Thus, the local 
AQC can only maintain its properties if W < g m . Since 
W ~ 1/Tdccoh, and tf ~ l/<7m for the local scheme in the 
weak-decoherence regime, the computation time is lim- 
ited by the decoherence </ < Td CC oh in the same way as in 
gate model QC. Therefore, the advantageous scaling of 
the local AQC requires phase coherence throughout the 
evolution as in the gate model. Insensitivity of AQC to 
decoherence only holds for the global scheme and does 
not apply to local AQC. It should be emphasized that in 
our treatment we have assumed that the minimum gap 
is a result of a first order quantum phase transition for 
which two-state model holds and the broadening of the 
energy levels and also thermal excitation do not mix the 
lowest two states with other excited states. For stronger 
noise or higher temperatures, one needs to take higher 
states into consideration. 
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W.M. Kaminsky, T. Lanting, D.A. Lidar, and G. Rose 
for fruitful discussions. 
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